Search for: All records

Abstract High-velocity clouds (HVCs) may fuel future star formation in the Milky Way, but they must first survive their passage through the hot halo. While recent work has improved our understanding of the survival criterion for cloud-wind interactions, few observational comparisons exist that test this criterion. We therefore present an initial comparison of simulations with the Smith Cloud (SC; d = 12.4 kpc, l, b = 40○, −13○) as mapped with the GALFA-HI survey. We use the Smith Cloud’s observed properties to motivate simulations of comparable clouds in wind tunnel simulations with Enzo-E, an MHD code. For both observations and simulations, we generate moment maps, characterize turbulence through a projected first-order velocity structure function (VSF), and do the same for HI column density with a normalized autocovariance function. We explore how initial cloud conditions (such as radius, metallicity, thermal pressure, viewing angle, and distance) affect these statistics, demonstrating that the small-scale VSF is sensitive to cloud turbulence while large scales depend on cloud bulk velocity and viewing angle. We find that some simulations reproduce key observational features (particularly the correlation between column density and velocity dispersion) but none match all observational probes at the same time (the large scales of the column density autocovariance is particularly challenging). We find that the simulated cloud (cloud C) showing growth via a turbulent radiative mixing layer (TRML) is the best match, implying the importance of TRML-mediated cooling for Milky Way HVCs. We conclude by suggesting improvements for simulations to better match observed HVCs. 

ABSTRACT The early growth of black holes (BHs) in atomic-cooling haloes is likely influenced by feedback on the surrounding gas. While the effects of radiative feedback are well-documented, mechanical feedback, particularly from active galactic nucleus (AGN) jets, has been comparatively less explored. Building on our previous work that examined the growth of a 100 $${\rm M_\odot }$$ BH in a constant density environment regulated by AGN jets, we expand the initial BH mass range from 1 to $$10^4\, {\rm M_\odot }$$ and adopt a more realistic density profile for atomic-cooling haloes. We reaffirm the validity of our analytic models for jet cocoon propagation and feedback regulation. We identify several critical radii – namely, the terminal radius of jet cocoon propagation, the isotropization radius of the jet cocoon, and the core radius of the atomic-cooling halo – that are crucial in determining BH growth given specific gas properties and jet feedback parameters. In a significant portion of the parameter space, our findings show that jet feedback substantially disrupts the halo’s core during the initial feedback episode, preventing BH growth beyond $$10^4 \, {\rm M_\odot }$$. Conversely, conditions characterized by low jet velocities and high gas densities enable sustained BH growth over extended periods. We provide a prediction for the BH mass growth as a function of time and feedback parameters. We found that, to form a supermassive BH ($$\gt 10^6 \, {\rm M_\odot }$$) within 1 Gyr entirely by accreting gas from an atomic-cooling halo, the jet energy feedback efficiency must be $$\lesssim 10^{-4} \dot{M}_{\rm BH} c^2$$ even if the seed BH mass is $$10^4 \, {\rm M_\odot }$$. 

Abstract The first infall of the LMC into the Milky Way (MW) represents a large and recent disruption to the MW circumgalactic medium (CGM). In this work, we use idealized, hydrodynamical simulations of an MW-like CGM embedded in a dark matter halo with an infalling LMC-like satellite initialized with its own CGM to understand how the encounter is shaping the global physical and kinematic properties of the MW CGM. First, we find that the LMC drives order-unity enhancements in MW CGM density, temperature, and pressure due to a $\mathcal{M}\approx 2$shock from the supersonic CGM–CGM collision. The resulting shock front extends from the LMC to beyond ∼R_200,MW, amplifying column densities, X-ray brightness, thermal Sunyaev–Zeldovich distortion, and potentially synchrotron emission from cosmic rays over large angular scales across the southern hemisphere. Second, the MW’s reflex motion relative to its outer halo induces a dipole in CGM radial velocities, withv_R ± 30–50 km s⁻¹atR > 50 kpc in the northern and southern hemispheres, respectively, consistent with measurements in the stellar halo. Finally, ram pressure strips most of the LMC’s CGM, leaving ∼10⁸⁻⁹M_⊙warm ionized gas along the past orbit of the LMC, moving at high radial and/or tangential velocities ∼50–100 kpc from the MW. Massive satellites like the LMC leave their mark on the CGM structure of their host galaxies, and signatures of such interactions may be observable in key all-sky tracers of the MW CGM and those of other massive galaxies. 

Abstract We explore the effect of variations in the Population III initial mass function (IMF) and star-by-star feedback on early galaxy formation and evolution using the Aeossimulations. We compare simulations with two different Population III IMFs:M_char = 10M_⊙and $M_{\max }=100M_{\odot }$(Aeos10) andM_char = 20M_⊙and $M_{\max }=300M_{\odot }$(Aeos20). Aeos20 produces significantly more ionizing photons, ionizing 30% of the simulation volume byz ≈ 14, compared to 9% in Aeos10. This enhanced ionization suppresses galaxy formation on the smallest scales. Differences in Population III IMF also affect chemical enrichment. Aeos20 produces Population II stars with higher abundances, relative to iron, of light andα-elements, a stronger odd–even effect, and a higher frequency of carbon-enhanced metal-poor stars. The abundance scatter between different Population II galaxies dominates the differences due to Population III IMF, though, implying a need for a larger sample of Population II stars to interpret the impact of Population III IMF on early chemical evolution. We also compare the Aeossimulations to traditional simulations that use single stellar population particles. We find that star-by-star modeling produces a steeper mass–metallicity relation due to less bursty feedback. These results highlight the strong influence of the Population III IMF on early galaxy formation and chemical evolution, emphasizing the need to account for IMF uncertainties in simulations and the importance of metal-poor Population II stellar chemical abundances when studying the first stars. 

Abstract The circumgalactic medium (CGM) of star-forming dwarf galaxies plays a key role in regulating the galactic baryonic cycle. We investigate how susceptible the CGM of dwarf satellite galaxies is to ram pressure stripping in Milky Way–like environments. In a suite of hydrodynamical wind tunnel simulations, we model an intermediate-mass dwarf satellite galaxy (M_*= 10^7.2M_⊙) with a multiphase interstellar medium (ISM;M_ISM= 10^7.9M_⊙) and CGM (M_CGM,vir= 10^8.5M_⊙) along two first-infall orbits to more than 500 Myr past pericenter of a Milky Way–like host. The spatial resolution is ∼79 pc in the star-forming ISM and 316−632 pc in the CGM. Our simulations show that the dwarf satellite CGM removal is fast and effective: more than 95% of the CGM mass is ram pressure stripped within a few hundred megayears, even under a weak ram pressure orbit where the ISM stripping is negligible. The conditions for CGM survival are consistent with the analytical halo gas stripping predictions in McCarthy et al. We also find that including the satellite CGM does not effectively shield its galaxy, and therefore the ISM stripping rate is unaffected. Our results imply that a dwarf galaxy CGM is unlikely to be detected in satellite galaxies; and that the star formation of gaseous dwarf satellites is likely devoid of replenishment from a CGM. 

Context.Feedback from stars in the form of radiation, stellar winds, and supernovae is crucial to regulating the star formation activity of galaxies. Dwarf galaxies are especially susceptible to these processes, making them an ideal test bed for studying the effects of stellar feedback in detail. Recent numerical models have aimed to resolve the interstellar medium (ISM) in dwarf galaxies with a very high resolution of several solar masses. However, when it comes to modeling the radiative feedback from stars, many models opt for simplified approaches instead of explicitly solving radiative transfer (RT) because of the computational complexity involved. Aims.We introduce the Realistic ISM modeling in Galaxy Evolution and Lifecycles (RIGEL) model, a novel framework to self-consistently model the effects of stellar feedback in the multiphase ISM of dwarf galaxies with explicit RT on a star-by-star basis. Methods.The RIGEL model integrates detailed implementations of feedback from individual massive stars into the state-of-the-art radiation-hydrodynamics code,AREPO-RT. It forms individual massive stars from the resolved multiphase ISM by sampling the initial mass function and tracks their evolution individually. The lifetimes, photon production rates, mass-loss rates, and wind velocities of these stars are determined by their initial masses and metallicities based on a library that incorporates a variety of stellar models. The RT equations are solved explicitly in seven spectral bins accounting for the infrared to He IIionizing bands, using a moment-base scheme with the M1 closure relation. The thermochemistry model tracks the nonequilibrium H, He chemistry as well as the equilibrium abundance of C I, C II, O I, O II, and CO in the irradiated ISM to capture the thermodynamics of all ISM phases, from cold molecular gas to hot ionized gas. Results.We evaluated the performance of the RIGEL model using 1 M_⊙resolution simulations of isolated dwarf galaxies. We found that the star formation rate (SFR) and interstellar radiation field (ISRF) show strong positive correlations with the metallicity of the galaxy. Photoionization and photoheating can reduce the SFR by an order of magnitude by removing the available cold, dense gas fuel for star formation. The presence of ISRF also significantly changes the thermal structure of the ISM. Radiative feedback occurs immediately after the birth of massive stars and rapidly disperses the molecular clouds within 1 Myr. As a consequence, radiative feedback reduces the age spread of star clusters to less than 2 Myr, prohibits the formation of massive star clusters, and shapes the cluster initial mass function to a steep power-law form with a slope of ∼ − 2. The mass-loading factor (measured atz = 1 kpc) of the fiducial galaxy has a median ofη_M ∼ 50, while turning off radiative feedback reduces this factor by an order of magnitude. Conclusions.We demonstrate that RIGEL effectively captures the nonlinear coupling of early radiative feedback and supernova feedback in the multiphase ISM of dwarf galaxies. This novel framework enables the utilization of a comprehensive stellar feedback and ISM model in cosmological simulations of dwarf galaxies and various galactic environments spanning a wide dynamic range in both space and time. 

Abstract We study the effects of active galactic nuclei (AGN) feedback on the Lyαforest 1D flux power spectrum (P1D). Using theSimbacosmological-hydrodynamic simulations, we examine the impact that adding different AGN feedback modes has on the predicted P1D. We find that, forSimba, the impact of AGN feedback is most dramatic at lower redshifts (z < 1) and that AGN jet feedback plays the most significant role in altering the P1D. The effects of AGN feedback can be seen across a large range of wavenumbers (1.5 × 10⁻³ < k < 10⁻¹s km⁻¹) changing the ionization state of hydrogen in the IGM through heating. AGN feedback can also alter the thermal evolution of the IGM and thermally broaden individual Lyαabsorbers. For theSimbamodel, these effects become observable atz ≲ 1.0. At higher redshifts (z > 2.0), AGN feedback has a 2% effect on the P1D fork < 5 × 10⁻²s km⁻¹and an 8% effect fork > 5 × 10⁻²s km⁻¹. We show that the small-scale effect is reduced when normalizing the simulation to the observed mean flux. On large scales, the effect of AGN feedback appears via a change in the IGM temperature and is thus unlikely to bias cosmological parameters. The strong AGN jets in theSimbasimulation can reproduce thez > 2 Lyαforest. We stress that analyses comparing different AGN feedback models to future higher precision data will be necessary to determine the full extent of this effect. 

Abstract We investigate how stellar feedback from the first stars (Population III) distributes metals through the interstellar and intergalactic medium using the star-by-star cosmological hydrodynamics simulation, Aeos. We find that energy injected from the supernovae (SNe) of the first stars is enough to expel a majority of gas and injected metals beyond the virial radius of halos with massM_dm ≲ 10⁷M_⊙, regardless of the number of SNe. This prevents self-enrichment and results in a nonmonotonic increase in metallicity at early times. Most minihalos (M_dm ≳ 10⁵M_⊙) do not retain significant fractions of the yields produced within their virial radii until they have grown to halo masses ofM_dm ≳ 10⁷M_⊙. The loss of metals to regions well beyond the virial radius delays the onset of enriched star formation and extends the period that Population III star formation can persist. We also explore the contributions of different nucleosynthetic channels to 10 individual elements. On the timescale of the simulation (lowest redshiftz= 14.3), enrichment is dominated by core-collapse supernovae for all elements, but with a significant contribution from asymptotic giant branch winds to thes-process elements, which are normally thought to only be important at late times. In this work, we establish important mechanisms for early chemical enrichment, which allows us to apply Aeosin later epochs to trace the evolution of enrichment during the complete transition from Population III to Population II stars. 

Abstract The Aeosproject introduces a series of high-resolution cosmological simulations that model star-by-star chemical enrichment and galaxy formation in the early Universe, achieving 1 pc resolution. These simulations capture the complexities of galaxy evolution within the first ~300 Myr by modeling individual stars and their feedback processes. By incorporating chemical yields from individual stars, Aeosgenerates galaxies with diverse stellar chemical abundances, linking them to hierarchical galaxy formation and early nucleosynthetic events. These simulations underscore the importance of chemical abundance patterns in ancient stars as vital probes of early nucleosynthesis, star formation histories, and galaxy formation. We examine the metallicity floors of various elements resulting from Population III enrichment, providing best-fit values for eight different metals (e.g., [O/H] = −4.0) to guide simulations without Population III models. Additionally, we identify galaxies that begin star formation with Population II after external enrichment and investigate the frequency of carbon-enhanced metal-poor stars at varying metallicities. The Aeossimulations offer detailed insights into the relationship between star formation, feedback, and chemical enrichment. Future work will extend these simulations to later epochs to interpret the diverse stellar populations of the Milky Way and its satellites.